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1 Introduction 



The most straightfoward process to analyze in QCD is the inclusive electron-positron an- 
nihilation into hadrons. The leading order QCD calculation confirms the charge and color 
assignments of the quarks, and gives a rough measure of the heavy quark masses. Further, a 
comparison of the next-to-leading order (NLO) calculation with experiment can in principle 
provide a direct measurement of the running of the strong coupling constant. 

In the perturbative QCD calculations, the final states that are summed over consist of colored 
quarks and gluons. It is assumed that the subsequent hadronization of these states affects 
the cross section only by a factor whose difference from unity is suppressed by powers of 
m\/q^, where m\ is the squared hadronic mass scale and is the invariant squared mass 
of the incident e^e~ pair. So far, this assumption has not been testable because of the 
uncertainties in the experiments. 

In this paper we analyze the analogous process in 1+1 dimensional massless quantum electro- 
dynamics, QED2. The model exhibits both confinement and asymptotic freedom and yet it is 
exactly solvable ^ |^, Q . Thus we may compare analytically the perturbative calculations 
with the exact results and study the effects of hadronization. What we find is an agreement 
between the NLO perturbative cross section and the corresponding term in the asymptotic 
expansion of the exact result (both of order g'^/q^ where g is the QEDg coupling constant). 
However, at higher orders in g^ jq^^ we find in the asymptotic expansion additional terms 
that depend also on Xivg^jq^. We did not expect that perturbation theory could generate 
these terms and we instead associate them with the hadronic mass logarithms. 



2 Exact and Perturbative Cross Sections 



We calculate the analogue of e^e~ hadrons in the following extension of the massless 
Schwinger model 0]. The Schwinger model fermion, t/^, is Yukawa-coupled to a massless 
scalar field, 0, (our "photon" ).0 A massless fermion / (our "lepton") with zero QEDg charge 
is also Yukawa-coupled to and for simplicity we will take both fermions to have Yukawa 
charge e. The full Lagrangian is 

L = FV4 + ^ + g /A)ij + fi jdf + 1/2000 + e(V^^ + /7)0. (1) 

So g is the analogue of the QCD4 strong coupling constant and e is the analogue of the 
electromagnetic coupling. We will confine all of our calculations to the lowest order in e and 
compare the exact (all orders in g) result with the result of the perturbation expansion in 
g. The exact results may be calculated using the dual realization of the Schwinger model in 

All of our amplitudes will involve a single factor of the momentum space propagator of this scalar field 
so we do not need to worry about the nonexistence of massless scalar fields in 1 + 1 dimensions. Giving this 
scalar small mass only effects the results by an irrelevant factor. 
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terms of a free scalar, h, of mass nih = g/ \pn and the bosonization formula Q 

= c : cos(20F/i) : . (2) 

The process that we consider is then // annihilation 

/(^i) + /(A:^) - X, 

where where the momentum assignments of the fermions are given in the parantheses and 
the final states X consist of quanta of the scalar field h because the massless Schwinger 
Model is equivalent to the free scalar free theory. The cross section at leading order in e^, 
but all orders in ^ is |[ 

^=4^^(^^)- (3) 

where g = /ci + /c2 is the energy of the virtual scalar photon, and -R(g^) is defined by 

/2(g2) = (i?+(g^) + i?4g2))/2 (4) 

and 

R±iQ^) = J (i^a;exp(iga;) exp(±47rAm^(x)), (5) 

Here Am^(a;) is the Wightman function for the free massive scalar field h and is calculated 
to be 

/dp 2 
, exp{ip ■ x) = —Ko{mhV-x'^ + zea;°). (6) 



and is the spherical Bessel function. 

The parton model calculation of // annihilation is performed perturbatively in g where the 
final states that are summed over consist of quanta of ip and of A^^. The Infra- Red (IR) 
collinear divergences that afflict QCD4 also occur in QEDg, but in 1 + 1 dimensions as power 
singularities rather than the logarithmic singularities of QCD4. We choose to regulate these 
singularities by temporarily giving the A'^ field a mass, mg. Then, as in QCD4, the KLN 
theorem |0 ensures that the inclusive cross section, summed over possible ip and of A^^ final 
states, will be finite as — at each order in perturbation theory. At leading order in g 
the only process that contributes is 

f{ki) + f{h)-^HPi) + HP2), 
and the cross section is easily calculated to be 

^The prefactor c is a normal ordering dependent constant and it equals c = (^2^)^ /2 '^hen the normal 
ordering mass equals m/j. 
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As in QCD4, the NLO corrections come from the interference of the leading order process 
with its one loop corrections, virtual gluon corrections, dal, and gluon bremstrahlung from 
one of the outgoing quarks, dal. 

First consider the virtual corrections. They consists of virtual correction to the vertex and 
self energy corrections of quarks. One of the reasons that we chose to regulate the IR coUinear 
divergence by giving the gluon a mass is that in 1 + 1 dimensions the one loop fermion self 
energy graph vanishes. This is because the amplitude is proportional to 

Trb^li*- Jk)^,] = (2 - d)Tr[i>- ^] = 0, (8) 

where is the fermion momentum and is the loop gluon momentum. Note that this trace 
vanishes For d = 1 + 1 spacctime. So we only need to consider the vertex loop correction 
of a gluon exchange between the outgoing quarks. A direct calculation gives the invariant 
amphtude for this process to be 

M,-eg J ai, ^ ^^^^^^ _ ^^^^^^ _ . ^yj 

Using the formula 7^ /i Pj'^ — 2 p /i for 1 + 1 dimensions and standard Feymann integration 
aproach one can prove that 



9> 

= -2^Ti^/2E^2{ml + q") 1^ dx dy 



X 



{1 — X — y)m?g — xyq 
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= -27ri7zB^4;-ln/3, (10) 
where /? = m'^/q'^. 

The interference of this one-loop amplitude with the leading amplitude then gives the fol- 
lowing order g^/q^ cross section: 

Next consider gluon bremstrahlung 

The invariant amphtude for this process is 



Mb = -ie^g^\uipi] 



^ Pl+ P3 _ P2+ p3 ^ 

{pi+psy {p2+P3y 



vip2). (12) 



Where is the "gluon" momentum and e its polarization. It is a simple exercise in kinematics 
to derive the corresponding cross section 

g"^ I 1 

do^l = TT^^^-. dxi (13) 

ZTT^* q^ Trpl — xi 
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where Xi = 2pi ■ q/q^. 



The corresponding gluon-emmission differential cross section in QCD4 is finite (see, for ex- 
ample, 1^), and the coUinear divergences arise only upon the integration over phase space. 
Indeed, apart from the phase space points where the gluon is collinear with one of the quarks, 
the differential cross section describes the probability of a 3-jet event which is observable and 
so the corresponding cross section must be finite. In contrast, in ci = 1 + 1 the emitted gluon 
is forced to be collinear, there are no jets, and so the differential cross section is infinite for 
/5 = 0. This is fine because the only physical observable is the total cross section which, as 
we will see, is finite. 

The maximum value of xi is 1 — /?, So integrating the differential cross section gives the NLO 
gluon emission cross section to be 

1 (7^ 1 

The result for the total cross section including including all of the order g"^ /q^ perturbative 
corrections is thus 



a = a 



' + ^1 + ^1=^' (15) 



That is, the loop correction, al exactly cancels the bremstrahlung correction, a^. The 
cancellation between these cross sections of the terms that are singular in is not a surprise 
- it is understood to be a consequence of the KLN theorem 0. However, it is surprising 
that, unlike QCD4 there is no finite residue. To compare this NLO result with the exact 
result in equation (H), we will next expand the functions R±{q^) in powers of /q^- 



3 Asymptotic Expansion of the Exact Result 

The functions R±{q^) which describe all of the interesting dynamics are compactly given in 
terms of 2-dimensional Fourier transforms (FT) of the exponential of the Wightman function 
in equation (|). It turns out that the evaluation of the asymptotic expansion of the R±{q^) 
may be reduced to a straightfoward exercise involving the 1-dimensional FTs of generalized 
functions [§. We will describe briefly the steps involved in the expansion bescause the 
calculation is rather involved. We will then evaluate explicitly the first few terms of the 
expansion and demonstrate the occurrence of the advertised \ng'^ /q^ term at order g^ /q^. 

We utilize the fact that when g — > 00, the FT of a function is dominated by its singularities 
on the real axis, in this case, at the origin. So all we need to do is to pick out the power 
and logarithmic singularities of the function exp [±47rAm^(a;)], and organize them according 
to their power. All the analytic terms will drop out, as their contributions are exponentially 
small. 

The first step in carrying out the asymptotic expansion of R±{q'^) is to expand the Wightman 
function that appears in the argument of the exponential in equation (|5|). The expansion of 
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Ajnix) given in equation is (this is the standard expansion of spherical Bessel functions 

i) 

oo 1 2 

47rA„,(x) = E ^i^K(-^' + ^ex'tmk + 1) - ln[^(-x^ + ^ex')]). (16) 

This expansion appears in the argument of the exponential in equation (|^). When the 
exponential is expanded and the terms in the double expansion are grouped according to 
their power dependence on m|(— + iex^). All we have to evaluate are 2-dimensional FTs 
of terms of a power times a logrithmic function. 

Next, by using the equalities 

- a;^ + zex° = -2(x+ - ze)(x~ - ze), (17) 

and 

(fx = dx^dx~ , (18) 

where 

= ^{x^ ±x^). (19) 

The expansion of R±{q^) is reduced to taking products over the 1-dimensional FTs of power 
and logrithmic functions of x"*" and x~ . 

dx{x^ - ieY ln"'(x^ - ie) exp ix^p, (20) 

-J 

where n and m are integers. 

We have carefully kept track of the regularizations which are need to make these FTs well- 
defined. The techniques are discussed in detail in the book by Gel'fand and Shilov [|]. 



The FTs in ( PDj ) may all be generated by taking derivatives with respect to A of the following 
FT equality given in entries 25 and 26 in the Table of Fourier Transforms in : 

/■°° , / , . X A • 27rexp (=FiA7r/2) , , 

j dx{x±ieY^expixp = ^j^^ —P^- (21) 

The subscript =F for p on the right-hand-side (rhs) of (^) defines generalized functions: 

= \p\>^-^e{Tp). (22) 
and 9 is the usual step function. To generate the desired FTs we expand A in Laurent 



series the terms on the rhs of (pi]) , collect terms with the same power dependence, and take 
derivatives with respect to A. The result is a series of FT of the functions of powers of xiie 
times powers of the logarithm of x ± ie. For example, ^ 



/CO 27r ( 

dx{x — ieY ~ 6xp ixp = — < 5"(p) 
-CXD t I 



^(n+l)-| 



2n\( 



(23) 



^ Here 'ipi^ + 1) is the logarithmic derivative of the usual gamma function evaluated at n + 1, so ^^{x) 
Refer to the section on the psi function in ||]. 
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and 



dx{x — ie)"'ln^(x — ie) exp ixp = — 



27r 
2n\(-)" 
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ilj'{n + 1) + V"^!^ + 1) - ^('^ + 1) - Y^^^ 



, [ln(p+) - ^(n + 1) + z7r/2] . (24) 

The 5 function terms drop out when \p\ is large. 

The final result of these three steps described above - expanding the integrand in the 2- 
dimensional FT definition of R±{q^) in powers of m\{—x'^ + zex°), factorizing each term into 
products of 1-dimensional FTs over x"*" and x~, and evaluating the resulting FTs using the 
master result in equation ( ^1]) - is an asymptotic expansion of R±{q^) in powers of g"^ /q^- The 
expansion can be taken out to arbitrary order by applying these steps. We have evaluated 
the terms out to order /q^ and find that for 7^ and > 0: 



R^i<i') = 1 + + ^^(12 + 41n[^]) + o(^), (25) 



and 



R,{q-) = ^1(6.96 - 3.79 ln[^]) + o(^). (26) 



There are several points to make regarding the above expansion. First, the fact that the terms 
of order vanish in both of the expansions in equations (25) and (p6| ) implies that the 



NLO term in the expansion of the exact cross section for // annihilation, (|^), also vanishes. 
This is in agreement with the perturbative result, equation ([T5|). Second, for either g"*" < 
or g~ < the functions R±{q^) vanish. Also, the R±{q^) contain delta function contributions 
concentrated on the lightcone, g"*" = or g~ = 0. Basically, all the terms in the FT's that 
are not delta functions depend on g^, see equation (pi]). Futhermore this property of R± - 
that it is nonzero only in the foward lightcone - persists at all orders in the expansion and 
can be traced back to the boundary conditions on the Wightman function Am(a;^) that leads 
to the regularization — > — iex^ . The third point is that the calculational steps outlined 
above imply that the asymptotic expansions of the R± functions have the form (again in the 
foward lightcone): 

00 2 2 2 

R^{<f) = 1 + E(^)'^K + alH^] + ■■■ + a:-'\n-'[^^]) (27) 

n=2 19 g 



and 



00 2 2 2 

R-il") = E i^n^l + hi ln[^] + ■ ■ ■ + hl'^ In-^'l^]), (28) 



n=4 



g' r g 



where the a's and 6's are g- and ^f-independent numbers. These expansions are all that we 
will need to know about the R± functions in this and succeeding work. For example, a review 
of the cross sections given in |T0[ together with the above result that the g"^ / q^ terms vanish 



imply that the NLO contributions to the cross sections must vanish. We will use this result 
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in subsequent work |Tl|] to compute the NLO distribution and fragmentation functions for 
the model defined above in equation (|l]). 



The final point is that, as advertised in the abstract, the asymptotic expansion of the exact 
result for the cross section contains a |6"ln[|2] term which comes from FT of ln^(a;^). This 
term is a new kind of non-perturbative effect, it is much more stonger than the common 
exponetially small non-pertubative effect ( e.g., instantons in QCD4 ) which are dropped in 
this expansion. Maybe it associated with hadronization. Note that for large enough the 
log term will dominate the other term (that we presume is calculable in perturbation theory) 
at order g^/q^. 



4 Conclusion 



Simple dimension counting in QCD4 distinguishes the perturbative corrections according to 
their twist. Thus, both the higher twist corrections and, by assumption, the hadronization 



2 



corrections are suppressed by powers of 1/q . The dimension counting is different in QED 
where the coupling g has the dimensions of a mass. Here all corrections are suppressed 
by powers of g'^/q^ oc m?/q^. Nevertheless, the NLO perturbative result agrees with the 
asymptotic expansion at order g^jq^. Further, the presence of the {g"^ / q^Y \n g"^ / q^ term in 
the expansion of the exact result is a surprise from the point of view of perturbation theory 
and poses a challenge to calculate. 



In the preceding paper [jTO[ the only property of the exact result that was used was that 
R{<f) — ^ 1 as ^ ±00. In this paper, the complete asymptotic expansion of R{q^) is 
described and the first few terms are calculated explicitly out to order ((7^/ q^Y. In subsequent 



work |TT[ we will use this expansion to extend the analysis in [|10| to the next-to-leading 
order. That is, to analytically check for process independence of the parton distribution and 
fragmentation functions in the massless Schwinger model. 
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